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Abstract 



We analyze the spatial and temporal resolving power of two-photon inten- 
sity interferometry for the light emitting source in single bubble sonolumines- 
cence (SBSL). We show that bubble sizes between several 10 nm and 3 jim 
can be resolved by measuring the transverse correlation function, but that 
a direct determination of the flash duration via the longitudinal correlation 
function works only for SBSL pulses which are shorter than 0.1 ps. Larger 
pulse lengths can be determined indirectly from the intercept of the angular 
correlator at equal photon frequencies. The dynamics of the bubble is not 

accessible by two-photon interferometry. 
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I. INTRODUCTION 



In single bubble sonoluminescence (SBSL) [|IJ an air bubble in water is trapped in the ve- 
locity node of an acoustical pressure field. Under defined conditions this bubble periodically 
emits intense, broad-banded flashes of light, synchronous to the driving sound. Measure- 
ments of pulse widths have given estimated values from less than 50 ps to more than 
250 ps ||, and an analysis of the bubble surface's temporal variation indicates minimal radii 
of ~ 0.5 yum PJ. A direct measurement of the size of the light emitting region does not exist 
so far. 

The fundamental light generating mechanism is still unknown. Some models |||| at- 
tribute the light emission to electronic excitations or Bremsstrahlung processes initiated by 
spherically converging shock fronts, leading essentially to a black-body spectrum. While 
reproducing the measured spectra quite well they require extraordinarily high temperatures 
with at least partial ionization of the trapped gas. The peak of this thermal spectrum and its 
high-energy Boltzmann tail are postulated to be hidden below the absorption edge of water 
(A < 180 nm) where the photon spectrum cannot be measured. Alternative explanations, 
invoking collision-induced emission || or quantum vacuum radiation [[?], operate at more 
moderate temperatures and, in the latter case, do not exhibit an (invisible) high energy 
component while still reproducing the spectral shape in the measurable low energy window. 
Numerical simulations @] combining the gas dynamics inside a sphere with the oscillations 
of the bubble surface yield pulse widths of ~ lps combined with a radius of ~ 0.5 /zm for 
the emission region. These simulations indeed produce strong shock waves which are seen 
as causal for the light emission, but the predicted flash duration does not seem to coincide 
with the experimental data. 

Intensity interferometry, based on Bose-Einstein correlations (BEC) between identical 
bosons, has been discussed as a possibility to clarify the structure and dynamics of the light 
emitting region PH by directly measuring its size and lifetime. This technique, originally 
developed to determine the angular size of stars |L(J, has recently been considerably refined 
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in order to extract the spatial and temporal structure of the hot reaction zone created in 
high energy nuclear collisions JTT]]. It has proven to be a valuable tool also in the analysis 
of dynamical processes in extremely small and short-lived particle sources. Its application 
to SBSL in the form of two-photon intensity interferometry should be further facilitated by 
the comparatively large number of photons emitted per pulse and by the absence of many 
of the complications present in high energy applications, like final state Coulomb and strong 
interactions among the produced particles and secondary production of particles outside the 
source by decay of unstable resonances. 

If successful, SBSL interferometry may in fact turn out to be a unique testing ground for 
the method itself since here, contrary to high energy physics where the sources are too small 
and shortlived to be probed externally, alternative investigation techniques are available 
which should allow for various cross checks. 

In this paper we supplement the suggestion of |8|U with a general and largely model 
independent discussion of the method and its application to SBSL, providing a quantitative 
analysis of its spatial and temporal resolving power in the limited range of experimentally 
accessible wavelengths. We show that the transverse correlator is sensitive to bubble sizes 
in the physically interesting domain while direct resolution of the pulse length via the lon- 
gitudinal correlator is probably not possible with present technology. It can be determined 
indirectly, however, via the intercept of the transverse correlator at zero opening angle 
between the detectors. The dynamics of the bubble can, unfortunately, not be resolved 
interferometrically. 



II. TWO-PHOTON CORRELATION FUNCTION 

The correlation function for two photons with momenta k a and fc& is defined as 

where Pi(k) is the inclusive single-photon spectrum and P2{k a , kb) is the two-photon coin- 
cidence spectrum. All photon energies are on-shell, u; a b = \k a b\- In the following we assume 



that the photons are emitted completely incoherently and that the source is spherically 
symmetric. While the latter is strongly suggested by the extreme stability of the oscillating 
bubble, some phase coherence among the emitted photons can presently not be excluded. 
Ideally one would test this by measuring the (true) intercept of the two-photon correlator 
([!]) at q = k a — k b = 0. As we will show such a measurement is difficult, due to the 
limited frequency resolution of present photon detectors; in the long run one might there- 
fore contemplate a comparison of 2- and 3-photon interferometry data to settle the issue of 
partial coherence |i2[| . 

For chaotic (incoherent) sources the single-photon spectrum P\{k) and the correlator 
C(k a , k b ) can be expressed in terms of the single-photon Wigner phase-space density S(x; K) 
of the emitting source fi3HT5f : 



p x {k) = / d 4 xS(x;k,uj), (2) 



C(k a , k b ) = 1 ' 
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Jd 4 xS(x;K,E) 

Here K = {oj a + u b , k a + k b )/2 and q = {uj a — u^, k a — k b ). The second equation in (|J) is 
an approximation in that the single-photon spectra in the denominator have been evaluated 
at the average momentum K rather than at k a and k b , and in both the numerator and 
denominator the correct energy variables (u> a , u b , and K Q = {uj a + u;&)/2, respectively) 
have been replaced by the on-shell energy corresponding to K, E = \k a + k b \/2. This 
approximation makes the following discussion more transparent, but can be systematically 
corrected for [|1^] (see below). The factor | in front of the second term in the correlator 



(0) takes into account that only photons with equal helicity states are affected by Bose 
symmetrization. 

Since the measured photons are on-shell and thus the Fourier transform in (||]) is not 
invertible, the space-time structure of S(x; K) cannot be uniquely reconstructed. Still, 
valuable information on the space-time structure of the source can be extracted from the 



measured correlation function in terms of the second central space-time moments of S(x; K) 
lEyiSyni . In the context of SBSL this will be discussed next. 



III. GEOMETRICAL ANALYSIS 

Detailed investigations (for a recent overview see have shown that the essential 
features of the correlator ([3]) can be captured by replacing the ^-dependence of the emission 
function S(x; K) by a Gaussian with the same center and width. This is even more true for 
SBSL applications than in high energy particle physics since here resonance decay effects 
which can invalidate this Gaussian approximation are absent. Inserting a Gaussian ansatz 
for S(x; K) into (|j) yields a correlator which is Gaussian in the relative momentum q. Due 
to the spherical symmetry of the problem there is only one distinguished direction which is 
defined by the photon pair momentum K. We therefore use a Cartesian coordinate system 
in which K = (E, 0, 0) (i.e. all independence can be expressed through the energy E) 
and q = (?h,5_l,0). We also have q = uj a — oj^ = - + q\\ ~ Following the techniques 



developed in ||i6| , |T8|| it is then easily seen that the most general form of the correlator reads 

C(q, £)wl + j e~ q - KX B )-«i ((*ir ct ~) 2 >( £ ) , (4) 

where the angular brackets denote averages taken with the source function, 

i ft \w77n J d 4 x f(x) S(x; E) 

(/(x))(E)= S d*zS(z;E) ' (5) 
and tilde superscripts indicate center-corrected coordinates, = — (xi)(E). (Note that 
(x±)(E) = due to spherical symmetry.) Eq. (P tells us that by measuring, at fixed E, 
the correlator as a function of q± and q n , respectively, we can determine the spatial variance 
(x 2 ± ) and the mixed variance ((ct — x n ) 2 ), respectively, of the effective source of photons with 
energy E. In principle, for different E the effective source can have different such "sizes" or 
"HBT radii" . 

In relativistic heavy ion collisions the .fT-dependence (here: ^-dependence) of the space- 
time variances plays an important role as a signature for collective expansion of the emitting 



source [|nj . This is most easily seen in the context of a hydrodynamically expanding, locally 
thermalized source whose momentum dependence is dominated by a boosted Boltzmann 
distribution ~ exp[— K ■ u(x) /T(x)} where u(x) is the collective expansion 4- velocity profile. 
This factor generates correlations between the momentum K and the position x in the 
emitter which in turn cause a independence of the HBT radii. The strength of these 
correlations can be estimated by writing 

S(x; E) ~ s(K-u(r, t)/T) « s{E/T){\ + 0{v/c)) . (6) 

Since the expansion velocity v(r,t) of the bubble in SBSL is limited by the shock velocity 
in the compressed bubble gas and thus below about 3 x 10 4 m/s Q (i.e. v/c < 10~ 4 ), the 
nonrelativistic estimate (|6]) is reliable and the x-K correlations induced by the collective 
dynamics of the bubble are weak. This is different for pion interferometry in heavy ion 
collisions: there the collective velocities are of the order of the light velocity, causing strong 
x-K correlations and an appreciable i^-dependence of the correlator which can be used as 
a diagnostic tool [ IT | . In SBSL interferometry, on the other hand, the independence of the 



HBT radii resulting from the weak x-K correlations is so small that it can be neglected 
in the measurable i?-range (see below). This is unfortunate since it means that SBSL 
interferometry will not give any direct access to the collective dynamics of the bubble during 
light emission |1| . On the other hand, it simplifies the theoretical description because we can 
neglect the ^-dependence of the HBT radii and also drop the cross term in the longitudinal 
HBT radius: 



;(x l - c *) a )«(^) + c '(P>. (7) 

This is true because both the displacement (x n ) of the source center in .Redirection and the 



cross term (a^ t) are also generated by the collective expansion |Tj| and thus here expected 
to be about 4 orders of magnitude smaller than the diagonal terms (x 2 ) and (t 2 ). 

The corrections resulting from the approximation (^) can now be systematically included 
following the discussion in ||16|| . Writing up to second order in q 
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(^^ + i^_^_*>!), (8 ) 

one derives to quadratic accuracy 

J <fx S{x- K, K ) e iqx « e AR2 ^ J d 4 x S{x; K) e iqx 

with 

A ^) = (^ lnP ' (B) )- (9) 

The single particle distribution Pi(u a ) may be similarly approximated by 

For P\{uJb) one obtains the same expression with the opposite sign for the term linear in q^. 
For the denominator in (|3|) one thus finds to quadratic order 

P x {u a )P x {u b ) « e ^*l«l+^\ j#(E) (11) 

with 

AJ2f(S) = ~lnP 1 (S). (12) 

Hence the corrected Gaussian expression (^) for the correlator reads 

C(q ± ,q ]l ,E)nl + ie- R ^ q2 ^ q t, (13a) 
R± = (xl) + AR\ , (13b) 
R 2 = (x 2 ) + c 2 (t 2 ) + AR 2 . (13c) 

Note that the corrections AR\ and AR 2 are proportional to the slope and curvature of 
the logarithmic intensity spectrum, respectively, and are therefore directly accessible from 
single-photon measurements. Both AR± and AR n turn out to be at most several 10 nm. 

The experimental realization of the correlation measurement (as proposed by Trentalange 
and Pandey ||) consists of two photo- multipliers focussing on the sonoluminating bubble 



at a relative angle 0. The required momentum resolution is achieved by suitably chosen 
apertures and pre-detector band-pass filters. The signal detected in one multiplier during 
one flash, proportional to the incident number of photons, is correlated with the output of 
the second device during the same flash and sampled over a sufficient number of bubble 
oscillations for statistics. 

This experimental setup suggests the use of q = u a — Ub and instead of the variables 
q± and q^. They are related via 

q\ = (aE> + - 2g 2 ) tan 2 §0 , (14a) 

9? = ll + (fl8 - Jj) tan 2 ±0 . (14b) 

R± and i?y can thus be isolated by fixing the average photon energy E and scanning the 
correlator either as a function of the opening angle at equal photon energies (go = 0, "trans- 
verse correlator"), or as a function of the energy difference go a t zero opening angle (0 = 0, 
"longitudinal correlator" ) : 

C(g = 0, 0, E) « 1 + § exp (-R\ AE 2 tan 2 \<p) , (15a) 
C(g Q , = 0, E) « 1 + | exp (-i? 2 g 2 ) . (15b) 

IV. RESOLVING POWER 

We will now proceed towards a quantitative estimate of the resolving power of such 
correlation measurements. We begin by noting that our approximations break down if the 
first terms on the right hand sides of (|13|b,c) become smaller than the corrections from 
the second terms. This turns out not to be the limiting factor, though, because similar 
lower limits for the HBT radii result from the fact that the opacity of water prohibits 
measurements at wavelengths in the ultraviolet below 180 nm, and low yields make the 
measurement difficult in the infrared (w ^5 1.5 eV). Practical measurements are only possible 
in the "transparency window" 1.5 eV ^ u ^ 6 eV (210 nm < A ^ 830 nm). As we will see 
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this limits, at fixed E, the opening angle <fi and the energy difference go which means that the 
correlator fllTSD can only be measured over a restricted interval of the control variables. If the 
correlator does not fall off appreciably over the accessible range, the HBT radius parameter 
cannot be accurately determined. This gives lower limits for R± and R n of several 10 nm, 
i.e. of the same order of magnitude as the upper limits for the correction terms Ai?y, AR±. 

After these general remarks let us enter a more detailed discussion, beginning with R±. 
In Fig. m we have plotted the correlator (|15a ) at fixed E = 3 eV and go — as a function 



of the opening angle 0. One sees that for R± = 10 nm the correlator falls off only by about 
20% over the measurable range, rendering the determination of R± difficult. According to 
Fig. [l|, good measurements of R± are possible for 10 nm < R± < 3 jttm. For R± > 3 //m 
angular resolution becomes a problem: the correlator falls off so rapidly that opening angles 
between the two detectors and angular apertures of each detector below 1° are required to 
resolve the correlation function. This obviously cuts down on event statistics. However, 



R± > 3 /im implies a source with transverse size y (x^) ^ 3 /im, see Eq. ( |13b| ). Since it is 
known that at the point of light emission the source is smaller than this (values below 1 /im 
are quoted in |1|), angular resolution of the measurement does not appear to be a crucial 
limiting factor. 

One should note, however, that go — as indicated in Fig. [I] implies ideal energy resolu- 
tion of the photon detector. We will see shortly that the finite energy resolution in real life 
modifies significantly the optimistic picture suggested by Fig. [l]. 

Let us now turn to a discussion of Ru . It is easy to see that if the light emitting source has 



a radius below 1 /im (i.e. y (x 2 ) ^ 1 /xm) then the r.h.s. of ( |13cj) is dominated by the duration 



of the light flash 5r = y (t 2 ) as soon as 5t becomes larger than about 3 femtoseconds (which 
corresponds to c5r = 1 /im). Since typical SBSL pulse durations discussed in the literature 
||T|-[3|] are much longer we can for the following discussion neglect in ( pL3c| ) the geometric 
contribution as well as AR± and write R^ ~ c5t. 

Fig. shows that SBSL pulses which last longer than 1 ps can only be resolved if the 
photon detector has an energy resolution well below 1 meV (!). This implies a relative band 



width 5X/X ?S 10~ 4 . As we will see in a moment, it is not sufficient that the filter-to-filter 
distance between the two detectors is known with this accuracy; the band width of each 
filter individually must satisfy this constraint. 

Commercially available filters in the visible region around 400 nm have band widths 
8X ^ 1 nm, corresponding to 5\/\ > 25 x 10~ 4 . For 5\ = 10 nm the authors of |3| quoted 
coincidence rates of 200-300 counts/s at a bubble-detector distance of 200 mm. For smaller 
S A the coincidence rate drops essentially like (<5A) 2 . 

Too large values of 5 A/ A imply that the correlator in Fig. |2| is averaged over a range qo 
which is much larger than the region over which the correlator drops back to 1 . This implies 
not only that the longitudinal correlation function C(qo) cannot be resolved, but also that 
the transverse correlation function C(0), being averaged over a wide go- ran g e > is strongly 
diluted. As a consequence C(4>) will not intercept the vertical (0 = 0)-axis at the ideal value 
1 + | = but at a much lower value. The effective intercept will be the smaller the larger 
the band width of the photon detector. 

It is not difficult to calculate the effective intercept value as a function of the ratio 
between the filter band width and the flash duration. Let us assume filters with a Gaussian 
frequency profile 

One easily checks that 

fu,a,5u>M fu> b ,Su>M = Sk ,8u>i{u) /g 0i2( w( A ^) (17) 

where 5uj' = 5u/\^2, cu = (cui + u)2)/2, Au = uj\ — tU2, and K = (u a + Ut>)/2, q Q = uj a — Ub 
as before. Neglecting x-K correlations in the source as discussed above we can assume that 
S(x;K) factorizes, S(x;E) w A(x) ■ s(E), where X(x) is normalized, / d A xX(x) = 1. For 
the product of single particle spectra in the denominator of Eq. ([|) we thus obtain 

P x {k a )Px{k h ) = 



(CO 



2 (5c;) 



(16) 



10 



J dUi du 2 fu a ,5u>(vi) fu> b ,5u(u2) s(uj 2 ) = 

fK ,5w'{u) J d(Au) f q0: 26u>(Au) 

xs(u + §Au/) s(u- |Aw) , (18) 



while the numerator is similarly found to be 



P 2 (k a ,k b )-P 1 (k a )P 1 (k b 



-j ' du (s(u)f f Kot 8u>{u) J d 4 xd 4 yX(x)X(y) 
xe - Mx - y) . {ea - eb) f d ( AuJ ) fqo ^,(Acu) 



x e 



iAw[(x°-» )-i(a>-v).(e„+e i )] 



where e a ^ are unit vectors in direction of k a ^. Since the filter band width Sou is narrow, 
the single particle spectrum s(E) can be taken constant inside the filter gap. This so-called 
"smoothness approximation" allows to perform the integration over Au in fll8~l) : 



P 1 {k a ) P 1 (fc 6 ) ^ J duo {s{uj)) 2 f Ko M") ■ ( 2 °) 

To obtain the effective intercept we divide (0) by ( (Of) and set k a = k b = K, i.e. go = — 
and e a = e^. Then the ^-integration in Eq. ( |19"1) factorizes, and the first factor on the r.h.s. 



of Eq. (|19D cancels against (20). The second factor can be easily evaluated in the Gaussian 
approximation where we replace the space-time factor X(x) by a Gaussian with the same 
rms widths. We find 

C(q = 0) = l + l 1 



^1 + A{5uo) 2 ((x 2 > + c 2 (t 2 >) /{Ticf 



1 1 

= ! + --_ (21) 

2 + 4(^) 2 J R 2 /(/ic) 2 

The dependence of this effective intercept on R n (which in the limit considered here is 
approximately equal to the flash duration c5r) is plotted in Fig. |3| for a fixed filter band 
width 8X = 1 nm at an average photon energy E = 3 eV (corresponding to A = 413 nm). 
For this case 5u R^/(hc) = 1 corresponds to a flash duration 5t = 88 fs ~ 0.1 ps. One 
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sees that the effective intercept is unity for St <C 0.1 ps; this is the domain where the 
given filter band width allows to resolve the flash duration by measuring the longitudinal 
correlation function ( |15b| ). For St ^> 0.1 ps the effective intercept decreases linearly with 
the pulse length, C(0) — 1 ~ I /St. Thus, assuming completely incoherent photon emission 



and cSt ^> y (j?), St can be determined from the effective intercept of the transverse two- 
photon correlator (|1 5a|) even if the longitudinal correlator (|15b|) does not show any structure 
for <5-values outside the experimental resolution (filter band width) Sou. 

It is worth pointing out that according to this analysis it is not necessarily advisable to 
strive for increasingly better filter resolution 8u resp. SX. As discussed above, increasing the 
band width enhances the coincidence rate quadratically while the effective intercept value 
decreases only linearly. In leading order the effect of SX on the experimental error bar of 
St thus cancels. Of course, determining St via the effective intercept of the correlator is a 
somewhat roundabout procedure which depends in a crucial way on the assumed chaoticity 
of the source; eventually one would like to achieve a genuine lifetime determination by 
measuring the longitudinal correlator ( |15b| ) with appropriate frequency resolution. 



V. SUMMARY 

We have shown that two-photon intensity interferometry can be applied to study the size 
and lifetime of the light-emitting region in single-bubble sonoluminescence. A measurement 
of the transverse correlation function in the experimentally accessible frequency range pro- 
vides sensitivity to sizes between several 10 nm and a few fim for the active bubble region. 
Present technological limitations on the frequency resolution limit a direct measurement of 
the flash duration via the longitudinal correlation function to pulse lengths below 0.1 ps. 
We showed, however, that for chaotic emitters with longer flashes the pulse duration can be 
determined indirectly via the intercept of the transverse correlation function as a function 
of frequency resolution. The dynamics of the bubble during light emission is not accessible 
by two-photon interferometry, due to the much too small expansion velocities. 
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FIG. 1. The correlation function C(E = 3eV, go = 0, as a function of £ = 2£'tan^0 for 
various values of R±. The domain above £ « 10.4 eV is not accessible due to light absorption in 
water. The second abscissa gives the detector settings for some typical £ values. 
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FIG. 2. The correlation function C(E = 3eV, qo, <j> = 0) as a function of the energy difference 
go, for various values of R±. One sees that pulse lengths above 1 ps require an energy resolution 
well below 1 meV. 
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